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We explore energy transfer in a generic three-level system, which is coupled to three non¬ 
equilibrium baths. Built on the concept of quantum heat engine, our three-level model describes 
non-equilibrium quantum processes including light-harvesting energy transfer, nano-scale heat trans¬ 
fer, photo-induced isomerization, and photovoltaics in double quantum-dots. In the context of light¬ 
harvesting, the excitation energy is first pumped up by sunlight, then is transferred via two excited 
states which are coupled to a phonon bath, and finally decays to the ground state. The efficiency of 
this process is evaluated by steady state analysis via a polaron-transformed master equation; thus 
a wide range of the system-phonon coupling strength can be covered. We show that the coupling 
with the phonon bath not only modifies the steady state, resulting in population inversion, but also 
introduces a finite steady state coherence which optimizes the energy transfer flux and efficiency. In 
the strong coupling limit, the steady state coherence disappears and the efficiency approaches the 
heat engine limit given by Scovil and Schultz-Dubois in Phys. Rew. Lett. 2, 262 (1959). 


I. INTRODUCTION 

With the rapid developments in measurement and ma¬ 
nipulation of microscopic systems, quantum effects such 
as coherence and entanglement are often utilized to en¬ 
hance the performance of microscopic devices. Even in 
biological systems, both experiments [I| and theoretical 
models [2, reveal that the long-lived quantum coher¬ 
ence may play an important role in highly efficient energy 
and electron transfer processes. How biological systems, 
such as light-harvesting complex, preserve such long-lived 
coherence and how nature benefits from the coherence 
are two key questions that define the emerging field of 
quantum biology. 

Taking a three-level system as a generic theoretical 
model, many interesting mechanisms can be well demon¬ 
strated and understood. Recently, the sunlight-induced 
exciton coherence is studied in a V-configuration three- 
level model Si- An interesting idea is to consider the 
energy transfer process from the perspective of heat en¬ 
gine (6|. For example, the coherence introduced by an 
auxiliary energy level can enhance the heat engine power 
[3, Hi- The early work considering a three-level maser 
model as a Carnot engine was carried out by Scovil and 
Schulz-DuBois [3, [13 ) yielding the heat engine efficiency 
rjQ and its relation with the Carnot efficiency. Later pa¬ 
pers elaborately reexamined the dynamics of this model 
by the Lindblad master equation and showed that the 
thermodynamic efficiency rjo is achieved when the output 
light-field is strongly coupled with the three-level system 
[lll - [l3 - The quantum heat engine provides us a heuris¬ 
tic perspective to better understand the basic physical 
processes in energy transfer and presents useful insight 
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to enhance the efficiency and output power in small sys¬ 
tems [ 1 M 3 . 

In this paper, we study the polaron effects of a phonon 
bath on the energy transfer flux and efficiency in a generic 
three-level model. The canonical distribution of a ther¬ 
mal equilibrium system requires a negligible coupling be¬ 
tween the system and its environment. As the coupling 
strength grows, the st ead^y state of the system will no 
longer be canonical (l§4^ . This non-canonical state ac¬ 
tually introduces the steady state coherence into the sys¬ 
tem without refereeing to specific forms of light-matter 
interaction or designing exotic system configurations. 
The bath-induced coherent effect is investigated by the 
polaron-transformed Redfield equation (PTRE) Ull > 
which bridges both the weak and strong system-bath cou¬ 
pling regions. The difference between the steady state 
efficiency and strong coupling limit tjq depends strongly 
on the phonon-induced coherence. Taking into account 
of the behaviour of both the flux and efficiency, we are 
able to optimize coupling and temperature in designing 
optimal artificial energy transfer systems. 

In this paper, we first introduce the three-level model 
and its non-equilibrium environment in section II, and 
then formulate the PTRE in section III. In section IV, 
the polaron effects of phonon-bath on the energy transfer 
flux and efficiency are studied in detail. We summarize 
our results in the last section. 


II. THREE-LEVEL SYSTEM MODEL 
A. Model system 

We consider the energy transfer process in the three- 
level system illustrated in FiglT] The site energy of the 
ground state |0) is set to zero. The two excited energy 
levels |1) and |2) form a two-level system (TLS, in the 
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following the TLS is referred to the two excited states), 
with the corresponding site energy ei and € 2 - The transi¬ 
tion due to the dipole-dipole interaction is characterized 
by J. Then the three-level system is modeled by the 
Hamiltonian Hq as: 

Ho = Y. e.|*)(*l + ^(|l)( 2 | + | 2 )(l|). ( 1 ) 

i=l,2 

We are interested in the transfer process in the single 
excitation subspace: The three-level system is firstly ex¬ 
cited to state | 1 ) by a photon held, then the excitation is 
transferred to state |2) through J (mediated by phonon 
modes), and hnally the excitation decays to the ground 
state |0) via spontaneous radiation. The pumping and 
trapping processes are modeled by the interaction with 
the two independent photon baths, which are coupled 
separately with two transitions | 0 ) O | 1 ) and | 0 ) t-)- | 2 ). 
The Hamiltonian of the photon baths and their interac¬ 
tions with the three-level system are given by 

Hi—p^t — ^ ^ T 1^) (^1 T H.C.^ , (2) 

k 

where ujik is the eigen frequency of the bath mode de¬ 
scribed by the creation (annihilation) operator (ciik), 
and its coupling strength to the excited state is gik- We 
note that the rotating wave approximation is applied in 
the system-bath interaction term. A phonon bath with 
creation and annihilation operators and bk of the bath 
mode Wvfc is coupled to the TLS via diagonal interaction 
with the coupling strength of fk- Thus, the phonon part 
is described by 

H. = Y + (|1) (1| - |2) (2|) ^ [fkbl + H.c.) . 

k k 

( 3 ) 

This microscopic three-level system immersed in the 
non-equilibrium environment was studied as a quantum 
heat pump phenomenologically without considering the 
details of the system-bath coupling 0 - In the case 
that the phonon bath is replaced by a single driving 
mode strongly coupled to the system, the dynamic steady 
states have been solved and the efficiency is given by 
m = £ 2/61 00 . In reality, the three-level model can 
be realized in both nature and laboratory. Taking the 
energy transfer process in photosynthetic pigment for ex¬ 
ample [Figl^a)], different baths could arise from differ¬ 
ent sources: the pumping light field (such as the sun-light 
photons) is considered as a high temperature boson bath; 
the trapping bath is formed by the surrounding electro¬ 
magnetic environment which models the energy transfer 
to the reaction center; and the phonon bath with in¬ 
verse temperature /3v describes the phonon modes cou¬ 
pled with the excited states. In addition, such a three- 
level (or more intermediate energy levels) system can be 
used to describe photoisomerization [Figl^b)], nanoscale 
heat transfer [Fig®c)] or photovoltaic current in 
double quantum dots [Figl5](d)]. 



Figure 1: (color online) The system is modeled by a three- 
level system: its ground state |0) and the excited state |1) 
(12)) is coupled with the pumping (trapping) bath; the ex¬ 
cited states |1) and |2) are diagonal-coupled with the phonon 
bath; the internal transition strength between |1) and |2) is 
characterized by J. The energy fluxes J'p, and Jt describe 
the energy exchange rate of the system with the pumping, the 
phonon and the trapping baths, respectively. The flux into 
the system is defined as the positive direction. 


(a) light-harvesting energy 
. transfer 
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(b) photoisomerization 




Figure 2: Realistic examples which can be studied by the 
three-level model with different heat baths, (a) In the photo¬ 
synthesis process, the three-level system works as an antenna 
that captures the energy from sunlight and then transfers to 
the reaction center, (b) Three eigenstates manifolds in pho¬ 
toisomerization. The bright states are pumped by the light 
held, then the populations relax to the intermedium and prod¬ 
uct states in the phonon environment, (c) The heat transfer 
in nanoscale can also use the three-level system as a bridge 
connecting the high temperature and low temperature heat 
baths, (d) In the electron transport problem, electrons tun¬ 
nel through double quantum dots which can be described by a 
three-level system. The quantum dot connects with a source 
and a drain. 


In this paper, we focus on the effects of the phonon 
modes on energy flux and efficiency. Usually when 
the system-phonon bath coupling strength is not weak, 
the Bloch-Redfield equation approach cannot be applied. 
Therefore, we will introduce the polaron-transformed 
Redfield equation (PTRE) [2^ , which gives reliable 

results from the weak to strong coupling region, to study 
the bath-induced coherent effects of this quantum sys- 








3 


tern. 


III. POLARON-TRANSFORMED REDFIELD 
EQUATION (PTRE) 


B. Definitions of energy fiux and transfer efficiency 


We are interested in the energy transfer flux and ef¬ 
ficiency of the three-level system at its non-equilibrium 
steady state. The steady state solution can be obtained 
by the master equation formally written as 

^^ ^ = (-^0 -b £p -b £v + Ps{t), (4) 

which describes the dynamics of the reduced density ma¬ 
trix (RDM) ps of the three-level system. The Liouville 
operator Lq denotes the non-dissipative term, £p, £v 
and Ct denote the dissipation effects associated with the 
pumping, phonon coupling, and trapping, respectively. 

To quantitatively investigate the energy transfer pro¬ 
cess, we define the steady state energy fluxes by calcu¬ 
lating the energy change of the three-level system: 

E{oo) = Tr^^idgJIt^oo = X! Trs [A [Ps{oo)] ids] 

i=p,v,t 

= Jdp -b jTv -b Jt- (5) 


It can be shown that Trg [CqHs] = 0. The three energy 
fluxes J7i, i = p, V, t are defined with respect to their cor¬ 
responding dissipation operator Ci. These energy fluxes 
have clear physical meanings of the energy exchange rate 
with the pumping field, phonon environment, and trap¬ 
ping field, respectively. In this work, we are interested in 
the steady state, in Eq.® the fluxes are calculated with 
Ps{oo), which is obtained by solving Ps{t) = 0. Straight¬ 
forwardly, we define the energy transfer efficiency by 


Jt(oo) 

Jp{oo) 


( 6 ) 


which is the ratio between the output and the input en¬ 
ergy fluxes. 

Without losing generality, we assume the pumping 
(trapping) bath is weakly coupled with the system and 
can be described phenomenologically by the local Liou¬ 
ville operator of the Lindblad form: 


^^[Ps] = ^[{n. + l){20-ps0+-{0t0-,ps}) 

+n, {20tpsO- - {0-0+, Ps})], (7) 


where f = p, t refers to the two photon baths, and 
rii are the corresponding decay rate and average pho¬ 
ton number, and the system operators are defined as 
Op = |1) (0|, 0+ = |2)(0|. The system-phonon bath 
coupling will be treated more rigorously as we are inter¬ 
ested in how this coupling affects the energy transfer over 
a broad range. To achieve this goal, we apply the PTRE 
equation, which will be introduced in the following sec¬ 
tion. 


A. Secular-Markovian Redfield equation in the 
polaron frame 


The Redfield master equation is valid up to the second 
order perturbation of the system-bath interaction. In or¬ 
der to go beyond this weak coupling limit, polaron trans¬ 
formation is introduced to incorporate the high-order 
system-bath interaction into the dynamics of the system. 
Here we focus on the coupling strength between the sys¬ 
tem and phonon bath, and the polaron transformation 
is only related to the two excited states. Therefore, it 
is convenient to consider the dissipative dynamics of the 
TLS first, then the resulting Liouville operator describ¬ 
ing the TLS dissipative process can be incorporated into 
the three-level system dynamics. We employ the Pauli 
matrix ax = |1) (2| -b |2) (1| and az — |1) (1| — |2) (2|, and 
define the polaron transformation 

=Ho + Hb + V, (8) 


where He = Hq + H„ is the Hamiltonian of the TLS 
with the phonon bath, the collective bath operator is 


Efc (fkbi - fkh) /Wvfc, and 


~ e J 

(9) 

Ho = -^cTz -b —Kax, 

Hb = 

(10) 

J 

(11) 

V = [o’x (cos B — k) + ay sin B]. 


The transformed system-bath interaction is V, where 
e = ei — £ 2 ) the expectation value of the bath oper¬ 
ator K = Trf, [pb cos B] is subtracted as a renormaliza¬ 
tion factor, and pi, is the thermal state of phonon bath. 
The spectrum function is chosen to be super-Ohmic 
as J (w) = 47rX)fe|/fc|^<5(w-a;fe) = 
where lOc is the cut-off frequency and a is a dimension¬ 
less parameter characterizing the system-bath coupling 
which is proportional to A/wc (A is the reorganization 
energy). Therefore we can obtain 


K = exp 


f 


J (w) I ! ^ 
du! -^ I riv (UJ) 


exp - 


1 - 


TTW^ 

2 





, ( 12 ) 


where n„{ui) = [exp(/3vw) — 1]“^ and ipi (x) = 

(n + x) ^ is the trigamma function. 

Since the thermal average of V is zero, then V is of 
the order of bath fluctuations and is a reliable perturba¬ 
tion parameter. Based on this consideration, the Born- 
Markov approximation is applied to derive the PTRE for 
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TLS in the Schrodinger picture as: 


dpe 

dt 


—i 


Pe 


^ pi^PeTpTa 


^pa'^aPe'Tp ^'^pTpPeTa]- 


(13) 


Here, pe is the RDM of the TLS in the polaron frame, 
and we use a new set of Pauli matrix Tq with respect 
to the eigenstates of the Hamiltonian Hq = e+ |+) (+| + 

e-l-) (-|: 


Tz = |+)(+|-|-)(-|, (14) 

T+ = !+)(-!, r_ = |-)(+|. (15) 

The corresponding eigenvalues and eigenstates are de¬ 
fined by 

e± = + {Kjf , (16) 

1+) = cos^ |1)-f sin^ |2), (17) 

I-) = sin^|l)-cos^|2), (18) 



Figure 3: (color online) The steady state of TLS as a func¬ 
tion of the system-bath coupling strength a. The steady 
states of the PTRE follow the canonical distribution in the 
polaron transformed basis, which rotates with the coupling 
strength a. In the weak coupling limit, the system steady 
state is the canonical distribution in the eigen basis (black 
dash line); while in the strong coupling limit, the steady state 
is the canonical distribution in the localized basis (red dot- 
dash line). The inset shows the coherent term of the steady 
state, which is small in the polaron transformed basis. We 
choose the parameters in units of J: ei/J = 5, £ 2 /^ = 4.5, 
uJc/J = 5 and /3vJ = 1. 


with tan0 = kJ/c. The transition rates T^^ are related 
to the half-side Fourier transformation of the bath corre¬ 
lation functions 

t2 pea 

ra/3 = ^y^ df(ec(±i)C/3(0)), (19) 

with 

iz{t) = sin 0 (cos i? (t) — k) , (20) 

S,±{t) = —[cos0 (cosH (t) — k) =F *sini? (<)](21) 


which are written in the form of the average values of 
the Pauli operators. The transition matrix M and the 
constant term (7^ = {Cz,Cx,Cy) are 



/ 7z 

'~)zx 

0 \ 


M = 

Ixz 

lx 

A -|- ^xy 1 , 

(23) 


\lyz 

^ lyx 

Iv ) 


II 

{C,,Cx,Cy), 


(24) 


The PTRE was firstly introduced by Silbey and cowork¬ 
ers [H, [ 2 ^, and has been widely used in solving the 
strong system-bath coupling problems. Moreover, it will 
be shown in Sec. IV that the results given by PTRE are 
consistent with those given by the Redfield equation in 
the weak coupling limit and the Eermi’s golden rule (or 
Forster theory) in the strong coupling limit jU, [ 2 ^, |2g . 
Therefore, the PTRE smoothly connects the two limits, 
and provides a useful tool to study the intermediate cou¬ 
pling region where there are usually no reliable approxi¬ 
mation methods. 


B. Steady state of PTRE 


For convenience, we rewrite Eq. (I13|) in the form of the 
Bloch equation 

|(r(t)), = -M(r(t)), + C. (22) 

Hereof (t))^ = [(r^ (t))^ , (r,, (t))^ , (r^ (t)) J with (^)^ = 
Trs[pe (t) ■] are the elements of the density matrix pe (t), 


where the eigenenergy level spacing is A = e+ — e_ and 
the expressions of the matrix elements are given in Ap¬ 
pendix A. 

The time evolution of (r (t))g is straightforwardly given 
by 


(r(t))^ = e-^* (t(0)),-M-iC +M-^C, (25) 


with the steady state {t{oo))^ = M~^C, and in the fol¬ 
lowing we will neglect time argument c» when referring 
to the steady state for convenience. The population dif¬ 
ference (tz)^ varies with the coupling strength as shown 
in EiglU In the weak coupling limit, the TLS steady state 
distribution is canonical with respect to its eigen basis, 
i.e.. 


_ I - exp{f3^Ve‘^ + .P) 
“^0 ^ I -I- exp(/3v\/e^ + J "^)' 


(26) 


which is also the result ~ exp(—/3v77o) given by equi¬ 
librium thermodynamics. When the system-bath cou¬ 
pling gradually increases, the system distribution devi- 
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ates from and follows the Boltzmann distribution 


1 — exp(/3v-\/e^ + K^,P) 

1 + exp{/3v-\/e^ + K^-P) ’ 


(27) 


with respect to the eigenenergy in the polaron trans¬ 
formed basis |-|-) and |—). In the strong coupling limit, 
we have 


lim 


(T-^)e 


1 - exp(/3ve) 

1 -I- exp(/3ve) ’ 


(28) 


which is the Boltzmann distribution with respect to the 
site energies ei and £2 of the localized basis |1) and |2). 
The deviation from the canonical state Pg™ due to the 
strong system-bath coupling has been studied via the cu- 
mulant expansion method in polaron transformed ther¬ 
modynamic distribution and from the view point 

of energy shell deformation HlEllii. 


C. Dissipative dynamics of the three-level system 
in the local basis 

Via the PTRE we can obtain a rather accurate descrip¬ 
tion of the TLS over a wide range of system-bath coupling 
strength. For further discussion on the property of the 
entire three-level system with the other two weakly cou¬ 
pled photon baths, all the observable quantities should be 
manipulated in the same frame of reference. To achieve 
this goal, we transform back into the frame of reference in 
the local basis. The diagonal terms of the system RDM 
are easy to deal with, as commutes with the polaron- 
transformation, 

(ctz (t)) = Trs+t, [ptot (t) CTz] = Trs [pe (t) a^] 

= Trs+t, [ptot (t) CTz] = Trs [pe (t) a^] , (29) 

where ptot (t) is the total density matrix for both the 
TLS and its bath, pe{t) = Tnlptot (t)], and ptotit) = 
^-i<7^b/2 (t) jg polaron-transformed total 

density matrix. However, the polaron transformation op¬ 
erator and ax{(Jy) do not commute with each other, thus 
the off-diagonal terms cannot be obtained exactly. We 
can use the approximation ptot (t) ~ pe (t) 0 Pb to obtain 
meaningful expressions for {ax (t)) and {ay {t)). This ap¬ 
proximation is essentially the Born approximation, which 
has already been used in deriving the PTRE. Based on 
these arguments, we have 


Following from Eas. (l^(l5T1)(l5^(l551) . the equations of 
motion for the TLS can be written in the form of 
[Pe{t)]ij = Z)m„[^v](ij.mn)[/3e(t)]rrm, then the expressions 
for the Liouville operator are straightforwardly ob¬ 
tained. 

The equations of motion of the three-level system are 
derived based on Eq.(|3]). The Liouville operator £v with 
polaron effects has been obtained from the PTRE of the 
TLS. One thing should be noted is that in the TLS, the 
conservation of population gives [pe{t)]ii + [Pe{t)]22 — 1, 
while in the three-level system the conservation relation 
becomes pm{t) + Pii{t) + P 22 {t) = 1, where Pij{t) = 
{i \ Ps{t) \j). The effects of the pumping and trapping 
baths are described by the Lindblad operator Cp and 
Ct defined in Eq.([7l)- Therefore, the PTRE for the three- 
level system is given as 


/ Pii{t) - P22{t) \ 
Pllit) + P22{t) I 

dt 3?[pi2(i)] I 
V 3[pi2(t)] / 


= -M 


( Pllit) - P22{t) \ 
Pllit) + P22{t) 

3^[P12(0] I 
V S[pi2(t)] / 
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/ IpUp - 7tnt 
JpBp + 7tnt 
0 

V 0 


(34) 


The matrix M is shown in Appendix B. The equations 
for the off-diagonal terms pmit) and pmit) are decou¬ 
pled from Ea. (|M)) and not related with the energy flux 
and transfer efficiency; thus poi it) and po 2 it) will not be 
involved in the following discussion. 


IV. ENERGY TRANSFER FLUX AND 
EFFICIENCY 

A. Steady state flux 

The steady state of the three-level system can be easily 
obtained from Ea. dMl) . which incorporates the polaron 
effects of the phonon bath. Then the steady state energy 
fluxes defined in Eq.([S]) are straightforwardly given as 

J'~y 

Jp = (^ilp[npPm - inp + l)pii] -( 3 ®) 

= e 27 t [^tPoo ~ i'^t + 1)^22]- i-Pi"^) 


{cTx it)) = KTrs [pe {t) axj, 

{ay it)) = KTr^ [pe it) ay]. (30) 

According to Eas. (fT7)l(fT51) . the Bloch vector {ait)) = 
[(cx it )), {ay {t)) , {az (t))] defined in the local basis of the 
TLS can be expressed with the quantities calculated in 
the polaron frame as 

(o-z(i)) = cos9{T:,it))^ + sin9{Txit))^, (31) 
(o-x(i)) = Ksin6 '(Tz(f))g-KCOs6»(Tx(t))g, (32) 
{ay it)) = -K{Tyit))^. (33) 


where we denote the steady state elements of RDM by 
Pij = (i|/9;,(oo)|j) for brevity. Figl3]presents energy fluxes 
with respect to a. In the extreme case that the system 
bath coupling is switched of! (a = 0), there is no loss of 
excitation energy, which results in Jp = —Jt, suggesting 
the input energy flux from the pump completely flows 
into the trap through the three-level system (note that 
we chose the positive direction as that the flux flows into 
the system). When the coupling turns on, a portion of 
energy flux leaks into the phonon bath thus Jp > —Jt- 
Both the pumping and trapping energy fluxes reach their 
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Figure 4: (color online) The steady state pumping (red solid 
line) and trapping (blue dashed line) energy fluxes versus a. 
The minus sign in front of the trapping flux suggests the en¬ 
ergy flows into the trapping bath. Both fluxes show a maximal 
value in the weak coupling case and then quickly decreases to 
zero when a increases. The inset shows the strong coupling 
case. We use the same parameters of the two level-system as 
in Fig|31 and the other parameters are chosen as: PpJ = 0.02, 
PtJ = 1, and 7 p/J = 7t/J = 0.01. 


optimal values in the intermediate coupling region and 
decrease to zero when the coupling strength is strong. 

In the context of heat engine, the trapping energy flux 
Jt in our model corresponds to the output power and 
corresponds to the input power. Usually, the power of 
a heat engine is small when the efficiency is high. Par¬ 
ticularly, at the maximal efficiency, all the processes are 
required to be quasi-static and take infinite time, and 
thus the power will be zero. To balance the conflict be¬ 
tween the efficiency and power, much work has been done 
to study the efficiency at maximum power (?7l - l^ . In the 
following, we will calculate the energy transfer efficiency 
of our system and show its competitive relation with the 
trapping flux, in analogy to the efficiency and power in 
the heat engine. 


B. Steady state efficiency 


Before presenting the result of efficiency defined in 
Eq.®, we begin with the analysis of the limiting cases. 
The first term on the right side of Eos. (1551) (IMl) depends 
only on the populations of the three-level system, and 
the second term represents the contribution of the off- 
diagonal terms (coherence in the local basis). As we 
have showed in Sec Iff, the steady state coherence in 
the local bases pi 2 vanishes in the strong coupling limit, 
then the efficiency is completely determined by the pop- 
nlations. According to the steady state solntion of the 
second equation in Eg. (1341) . we obtain the relation 


- {n^ + l)pii] = + l)p22 - ntpoo]. (37) 



Figure 5: (color online) (a) The steady states efficiency p, (b) 
the excited states population pn — p 22 and coherence K[pi 2 ] 
versus the system-phonon bath coupling strength character¬ 
ized by dimensionless parameter a. The dashed line indi¬ 
cates the strong coupling limit po in (a). When the popu¬ 
lations are inverted, p is less then po, the red dots indicate 
the corresponding range of a. The results given by the Red- 
field equation and the Fermi’s golden rule are shown with the 
dashed-dot lines. The strong coupling regions are plotted in 
the insets. All the parameters are chosen as the same as in 

FigH 


With this relation, the efficiency in the strong coupling 
limit reads 

^ _ e27t[(nt + l)p22 - ntpoo] ^ ez 
ei^p[{np + l)pii - Uppoo] ei' 

This result indicates that when the coherence is negligi¬ 
ble due to the strong system-phonon coupling, the energy 
transfer efficiency p approaches po, which is consistent 
with the key result of Refs.fl^. We notice that Eg. (1571) 
shows that the net rate of pumping one excitation to |1) 
equals to the net rate of trapping one excitation from 
12) to |0). In general, the efficiency is closely related 
to the phonon bath induced coherence [s^ of the ex¬ 
cited states. If we require the system outputs positive 
energy, i.e., 7 t(nt + l)p 22 > "ttntPoo, then according to 
Eos. (1^(1551) (1571) . 5i[pi2] > 0 leads to p > po and vise 
versa. 

According to our discussion of the flux in the last sub¬ 
section, when the coupling strength a = 0, the energy 
transfer efficiency p = 1 because there is no loss of energy 
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flux. When the coupling strength gradually increases, 
the efficiency decreases. However, after reaching its min¬ 
imum value, the efficiency starts to rise with a, which 
is shown in FiglSJa). The increase of efficiency assisted 
by noise was studied extensively in the context of energy 
transfer in light-harvesting systems As we fur¬ 

ther increase a, the efficiency grows beyond the strong 
coupling limit rjQ and then gradually approaches this limit 
from above. The strong coupling region is plotted in the 
inset of FigElJa). 

Interestingly, we find population inversion of the two 
excited states in the regimes of 77 > 779. We plot the pop¬ 
ulation difference between states |1) and |2) in FiglS^b). 
In the intermediate coupling region indicated between 
the two red dots, the steady state population satisfies 
Pii < P 22 (the effective temperature associates with these 
two states is positive), the corresponding efficiency rj is 
less then rjo as shown in FiglSJa). On the contrary, out¬ 
side this intermediate region, i.e., when the coupling is 
either very weak or very strong, the populations are in¬ 
verted pii > P 22 (the effective temperature is negative); 
meanwhile rj increase beyond 779 . In the local basis, the 
population and coherence are coupled with each other 
due to the polaron effects: The population inversion 
happens when 5 i[/ 3 i 2 ] < 0 [FiglUb)]. The fact that the 
population and coherence in the local basis have similar 
behaviour can be explained from Eg. (1311) and Eg. (1321) . 
Here, the coherence {tx (t))g in the polaron basis is negli¬ 
gibly small (see the inset of Figl3]) to have significant ef¬ 
fects, then the terms (cr^) = pn — P 22 and {ax) = 25R[/9i2] 
are both determined by {tz) 

In FigEKa), we also compare the efficiency 77 calculated 
by the PTRE method with those predicted by the Red- 
field equation and the Fermi’s golden rule approaches. 
As we mentioned before, in the weak and strong cou¬ 
pling limits, the PTRE method agrees with the Redfield 
equation and the Fermi’s golden rule, respectively, and 
it connects these two limits with a non-trivial minimum 
which is related to the coherence in the local basis. 


C. Further discussions 

1. kinetic models 

In the strong coupling regime, we can map this energy 
transfer process into a simple excitation kinetic model as 
shown in FigjSKa). Each step of energy transfer is de¬ 
scribed by an effective flux (J7p®, and The ef¬ 

fective transfer flux between the two excited states is 
approximately proportional to 72 [Figl^b)], which char¬ 
acterizes the relaxation rate of the two excited states. 
When (or 7 ^) is smaller than the trapping flux 
(or 7 t), the excitation in excited states will be quickly 
captured by the trapping field without enough time to 
first get equilibrated with the phonon bath. Conse¬ 
quently, the populations of the two excited states are 
inverted and the real part of the coherence becomes neg- 
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Figure 6 : (color online) (a) The kinetic model of a single 
excitation transfer cycle. In the strong coupling regime, the 
energy transfer processes between different local states can be 
described by the effective fluxes, (b) The transfer rate 7 ^ ver¬ 
sus system-bath coupling strength a. The red line indicates 
the trapping rate 7 t. The parameters here are the same with 

FigH 

ative. This phenomenological mechanism explains why 
the efficiency 77 is higher than 779 in the strong coupling 
limit. 

When the system-bath coupling strength becomes 
weaker, the local basis frame is no longer a good option 
for the kinetic picture. The two excited states couple 
with each other and can be together considered as an 
excited state manifold, as shown in Figl7l(a). The sin¬ 
gle excitation carrying certain amount of energy passes 
through the excited states | 1 ) and | 2 ), and its average 
residence time (t) in the excited states is negatively corre¬ 
lated with the transfer effici ency (in analogy to the light¬ 
harvesting efficiency in Ref. (3a. |34|1: i.e., the longer the 
excitation stays in the excited states, the more energy 
will be lost to the phonon bath, and the lower energy 
transfer efficiency will be. During a cycle that the single 
excitation starts from | 0 ) and finally returns to | 0 ), the 
average residence time {t) is proportional to the excited 
states population pn + P 22 at the steady states, as shown 
in FiglT^b). Though not quantitively exact, this kinetic 
model qualitatively explains the local minimal of the ef¬ 
ficiency 77 via the average residence time (t) pn -I- P22- 

2. temperature dependence 

Besides the system-phonon bath coupling strength, the 
temperature of the phonon bath also affects the energy 
transfer process, as shown in the two-dimensional con¬ 
tours of energy transfer efficiency [FiglSKa)] and trapping 
energy flux [Figl^b)]. The efficiency behaves the same 
at the high phonon bath temperature as in the strong 
coupling. In the high temperature limit, even when the 
coupling strength is weak, the efficiency is still close to 
779 . As seen from Eg. (1121) . in either limit a —> 00 or 
/3v —^ 0 , the renormalization factor k —^ 0 ; therefore, ex- 
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V. CONCLUSION 



Figure 7: (a) The two excited states form a black box for 
the input and output excitation due to the internal coupling, 
(b) The average residence time (t) is proportional to the total 
population of the excited states. The parameters here are also 
the same as those in Fig|4] 
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Figure 8: (color online) The dependence of (a) the energy 
transfer efficiency and (b) the trapping energy flux on the 
coupling strength and temperature of the phonon bath. The 
temperatures and dissipation coefficients of the pumping and 
trapping bath are the same as in the Fig[5] 


In this paper we use the polaron transformed Redfield 
equation (PTRE) to analyze the effects of the phonon 
bath on the energy transfer process in a generic three- 
level model. As a quantitative method, the PTRE can 
reliably describe the dependence of the steady state co¬ 
herence on the system-bath coupling strength ranging 
from the weak to strong coupling regime. Our analysis 
shows that the steady state coherence between the two 
excited states is crucial to the energy transfer efficiency. 
When the effective temperature of the excited states is 
negative (populations are inverted), the coherence carries 
a positive real part and enhances the efficiency beyond 
the strong coupling limit rjo- On the contrary, if the 
effective temperature is positive (populations are not in¬ 
verted) , the coherence carries a negative real part and is 
detrimental to the efficiency. The energy flux and effi¬ 
ciency compete with each other and cannot reach max¬ 
imum simultaneously; however, the study of their be¬ 
haviours with respect to the coupling strength and tem¬ 
perature provides the key information about how to make 
an optimal compromise between the two quantities. We 
will consider how to use quantum control to optimize the 
energy transfer process in the future study. 
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cept for the weak coupling and low temperature case, the 
efficiency 77 does not change obviously. 

The trapping energy flux has a different temperature 
dependences for weak and strong system-bath couplings. 
The flux —Jt grows (goes down) with descending /3v 
when a is small (large). Moreover, does not sen¬ 
sitively depend on with small a in contrast with the 
efficiency. When the coupling a is around 1, the flux — J7t 
changes no more than 20 % in amplitude comparing with 
its maximum. The optimization of the efficiency and the 
trapping flux can be achieved in two different regimes: 
1) The coupling strength is weak and the temperature 
of the phonon bath is high. 2) The coupling strength 
is medium (a ^ 2.5) and the bath temperature is low 
(/3v>l)- The first regime corresponds to the high temper¬ 
ature classical limit, and the second regime corresponds 
to low-temperature quantum regime, where bath-induced 
coherence enhances the energy transfer process. 


Appendix A M and C in Eg. (1231) and Eg. (1241) 

The quantities defined in Ea. (|23)) and Ea. (|24)) are de¬ 
termined by the superposition of the correlation functions 
Ea. dTHl) following Fas. (11^1) (1221) . The straightforwardly 
calculation gives: 


II 

M 


II 


Ty = 

2 

"fzx = 

4 

'^xz — 

4 


1 

sin20 J dtf {t) cos (At). 
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Ixy 


lyx 


Ivz 

C, 


a 

Cy 


where 


1 

/ dtg {t) sin {At), 

2 Jo 

—-K^J^ cos^ 9 / dtf {t) sin (At). 

2 Jo 

-K^J^sin20 / dtf (t) sin (At), 

4 Jo 

i r°° 

/ dtsin{At) 

^ J — CO 


X [cos^ 6 cosh [Q (t)] + sinh [Q 

/ OO 

dt sin {At) cosh [Q (t)], 

-OO 

—-K^J^sin20 / [1 — cos (Ai)] 

X [cosh [Q (i)] — cosh [Q (—i)]], 


/ {t) = cosh [Q (t)] + cosh [Q (—i)] — 2, 
g (t) = sinh [Q (t)] + sinh [Q (-i)] . 

Using the super-Ohmic spectrum J (w) 
the function Q {t) reads 


Q{t) = 


f 

JQ 


duj 


J (w) 


TTW^ 


[(2nv (w) + 1) cos (wt) — i sinwt] 


— 1 + ujit^ — i2uict 
= a ^-_ _ r, -h 


25R 


^1 ( 


1 

/HvCJc 




(/3vWc) 


Appendix B M in Eg. (1341) 

The Liouville operator for the three-level system is 
obtained from Ea. (E^ for the TLS with the expressions 


in Appendix A. Here the relation poo + Pii + P22 = 1 
should be used to substitute [pe]ii + [Pe]22 = 1 for the 
TLS. Taking the contributions of the Lindblad terms £p 
and Ct defined in Eq. ([ 7 ]) into consideration, the elements 
of the matrix M in Eq. (IMl) are 

Mil = Iz cos^ 6» -k 7a; sin^ 0 -k ^ (7a;z -k ^zx) sin 20 

+ 1 ) -k7t(nt + 1 )], 

Mi2 = -Cz COS0 - Ca;Sin0-I-i[7p(3np + 1) - 7t(3nt-I-1)], 

Mi3 = K~^Yixz sin^ 0 - ^zx COS^ 0 + ^(7^ - lx) sin20], 

Mi 4 = -k“^(A-|- 7a;y)sin0, 

M21 = 2 [ 7 p(^^p + 1 ) - lt{nt + 1 )], 

M 22 = ^[7p(3?rp + 1) + 7t(3nt-I-1)], 

M 23 = M 24 = 0 , ^ 

^31 = K,[^zx Sin^ 0 - -fxz COS^ 0 + 2 (7z - 7a:) siu 20], 

Mo 2 = k{Cx COS 9 — Cz sin 0 ), 

Moo = lx cos^ 0 -f 7z sin^ 0 - i(73;^ -k -jzx) sin 20 

+ ^[7p(np + 1) +7t(nt + 1)], 

M34 = (A -k 7a;y) COS 0 , 

M41 = K[(A-7ya;)sm0-7y2COS0], 

M42 — KCy , 

M43 = - (A - 7ya;) COS 0 - 7j,z sin 0 , 

M44 = ly + ^[ 7 p(np + 1) + 7 t(nt + !)]■ 
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